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XMM-Newton observation of the most X-ray-luminous galaxy 
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Abstract. We report on an XMM-Newton observation of RX J1347.5-1 145 (z=0.451), the most luminous X-ray cluster of 
galaxies currently known, with a luminosity L x = 6.0 + 0.1 x 10 45 erg/s in the [2-10] keV energy band. We present the first 
temperature map of this cluster, which shows a complex structure. It identifies the cool core and a hot region at radii 50-200 
kpc to south-east of the main X-ray peak, at a position consistent with the subclump seen in the X-ray image. This structure 
is probably an indication of a submerger event. Excluding the data of the south-east quadrant, the cluster appears relatively 
relaxed and we estimate a total mass within 1.7 Mpc of 2.0±0.4 x 10' 5 Mq. We find that the overall temperature of the cluster 
is kT = 10.0 ± 0.3 keV. The temperature profile shows a decline in the outer regions and a drop in the centre, indicating the 
presence of a cooling core which can be modelled by a cooling flow model with a minimum temperature ~2 keV and a very high 
mass accretion rate, M ~ 1900 MQ/yr. We compare our results with previous observations from ROSAT, ASCA and Chandra. 
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1. Introduction 

In this paper we present the first results from an XMM-Newton 
observation of RX J1347.5-1 145, the most X-ray-luminous 
galaxy cluster known (Schindler et al. 1995). This cluster has 
been detected in the ROSAT All-Sky Survey and further stud- 
ied with ROSAT HRI and ASCA (Schindler et al. 1995, 1997). 
It shows a very peaked X-ray emission profile and presents a 
strong cooling flow in its central region. Submm observations 
in its direction showed a very deep SZ decrement (Komatsu 
et al. 1999, 2001; Pointecouteau et al. 1999, 2001). Due to the 
presence of gravitational arcs, this cluster is also well suited for 
a comparison of lensing mass and X-ray mass. Optical stud- 
ies of weak lensing have been performed by Fischer & Tyson 
(1997) and Sahu et al. (1998). Recent Chandra observations 
(Allen et al. 2002) discovered a region of relatively hot, bright 
X-ray emission, located approximately 20 arsec to the south- 
east of the main X-ray peak at a position consistent with the 
region of enhanced SZ effect. This could be an indication for a 
subcluster merger, pointing to a complex dynamical evolution. 
A comparison of the XMM-Newton and SZ results, a more de- 
tailed analysis of the complex dynamical state of the cluster and 
a comparison of lensing mass and X-ray mass will be presented 
in a forthcoming paper (Gitti et al. in prep.). 

RX J1347.5-1145 is at a redshift of 0.451. With H = 
70kms~ 1 Mpc -1 , and fl^ = 1 - Qa = 0.3, the luminosity 
distance is 2506 Mpc and 1 arcsec corresponds to 5.77 kpc. 



Send offprint requests 
myriam . gitti@uibk . ac . at 



to: 



Myriam Gitti, e-mail: 



2. Observation and data preparation 

RX J1347.5-1 145 was observed by XMM-Newton in July 2002 
during rev. 484 with the MOS and pn detectors in Full Frame 
Mode with THIN filter. We used the SASv6.0.0 processing 
tasks emchain and epchain to generate calibrated event files 
from raw data. Standard processing was applied to prepare data 
and reject the soft proton flares. The remaining exposure times 
after cleaning are 32.2 ks for MOS1, 32.5 ks for MOS2 and 
27.9 ks for pn. The background estimates were obtained us- 
ing a blank-sky observation consisting of several high-latitude 
pointings with sources removed (Lumb et al. 2002). The back- 
ground subtraction was performed as described in full detail in 
Arnaud et al. (2002). The source and background events were 
corrected for vignetting using the weighted method described 
in Arnaud et al. (2001). 

3. Morphological analysis 

The adaptively smoothed, exposure corrected image 
(MOS+pn) in the [0.9-10] keV energy band is presented 
in Fig. ^ The image is obtained from the mosaic of the raw 
images corrected for the mosaic of the exposure maps by 
running the task asmooth set to a desired signal-to-noise 
ratio of 20. A number of notable features are visible. In 
particular, we note a sharp central surface brightness peak at 
13 h 47 m 30 s .6 - 11°45'09".0 (J2000), in very good agreement 
with the optical centroid for the dominant cluster galaxy 
(Schindler et al. 1995). We also confirm the presence of a 
region of enhanced emission ~ 20 arcsec to the south-east 
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Fig. 1. Total (MOS+pn) XMM-Newton EPIC mosaic image of 
RX J1347.5-1145 in the [0.9-10] keV energy band. The im- 
age is corrected for vignetting and exposure and is adaptively 
smoothed. 



(SE) of the X-ray peak and, on large scale (~ 80 arcsec), the 
extension of the X-ray emission to the south, already revealed 
in previous observations with Chandra (Allen et al. 2002). In 
Fig. |2 we show an overlay of the VLT image of the central 
region of RX J 1347.5 -1145 with the X-ray contours derived 
from Fig.^ 




radius (arcmin) 

Fig. 3. Background subtracted, azimuthally averaged radial 
surface brightness profile for SE quadrant data in the [0.3-2] 
keV range. The dotted line shows the profile in other directions 
(undisturbed cluster), which appears relatively regular and re- 
laxed. An excess of emission in the SE quadrant between radii 
of ~ 20-50 arcsec (100-300 kpc) is visible. 



The surface brightness profile of the undisturbed cluster 
is fitted in the CIAO tool Sherpa with various models, which 
are convolved with the XMM-Newton PSF. A single /?-model 
(Cavaliere & Fusco Femiano 1976) is not a good description 
of the entire profile: a fit to the outer regions (350 kpc - 1.7 
Mpc) shows a strong excess in the centre when compared to the 
model. The peaked emission is a strong indication for a cooling 
core in this cluster. We found that for 350 kpc- 1.7 Mpc the data 
can be described by a y6-model with a core radius r c = 367 ± 3 
kpc and a slope parameter/? = 0.93+0.01, while for r <350 kpc 
the data can be approximated by a jfi-model with r c — 40 + 0.2 
kpc and p = 0.55 + 0.02 (90% confidence levels). 





Fig. 2. VLT image of the central region of RX J 1347.5-1 145 
(Erben et al. in prep.). Superposed are the ([0.9-10] keV) XMM 
X-ray contours (levels: 0.003, 0.015, 0.045, 0.06, 0.15, 0.3, 0.6, 
1.5, 3 cts/s/arcmin 2 ). The image is ~ 4.4 x 4.3 arcmin 2 (North 
is up, East is left). 

We compute a background-subtracted vignetting-corrected 
radial surface brightness profile in the [0.3-2] keV energy band 
for each camera separately. The profiles for the three detectors 
are then added into a single profile, binned such that at least 
a sigma-to-noise ratio of 3 was reached. The cluster emission 
is detected up to 1.7 Mpc (~ 5')- In Fig. [3] we show the X-ray 
surface brightness profiles for the disturbed SE quadrant com- 
pared to that from data excluding the SE quadrant. We note that 
the data excluding the SE quadrant (hereafter undisturbed clus- 
ter) appear regular, while those for the SE quadrant (containing 
the X-ray subclump) show a clear excess of emission between 
radii of ~ 100 and 300 kpc relative to other directions. 



4. Temperature map 

The temperature image of the central cluster region shown in 
Fig.0]is build from X-ray colours. Specifically, we produce the 
mosaics of MOS images in four different energy bands ([0.3- 
1] keV, [1-2] keV, [2-4.5] keV and [4.5-8] keV), subtract the 
background and divide the resulting images by the exposure 
maps. A temperature is obtained by fitting the values in each 
pixel with a thermal plasma. In particular we note that the very 
central region appears cooler than the surrounding medium and 
the SE quadrant, which corresponds to the subclump seen in 
the X-ray image, is significantly hotter than the gas in other 
directions. 



5. Spectral analysis 

For the spectral analysis we treat the SE quadrant containing 
the X-ray subclump separately from the rest of the cluster. The 
data for the undisturbed cluster are divided into the annular re- 
gions detailed in Table 1. A single spectrum is extracted for 
each region and then regrouped to contain a minimum of 25 
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9.00 < — Temp (keV) --> 17.40 

Table 1. The results from the spectral fitting in concentric annular 
regions (undisturbed cluster). Temperatures (kT) are in keV, metallic- 
ities (Z) in solar units and [2-10] keV luminosities (L x ) in units of 
10 44 erg s . The total x 1 values and numbers of degrees of freedom 
(DOF) in the fits are listed in column 5. Errors are 90% confidence 
levels (Ax 2 = 2.71) on a single parameter of interest. 




Right Ascension (2000) 

Fig. 4. Temperature map obtained by using 4 X-ray colours 
([0.3-1], [1-2], [2-4.5], [4.5-8] keV) and estimating the ex- 
pected count rate with XSPEC for a thermal MEKAL model, 
with fixed Galactic absorption Nh = 4.85 x 10 20 cirT 2 and 
metallicity Z = 0.3Zq- Superposed are the X-ray contours. The 
features outside the last contours are not significant, as they are 
mainly due to noise fluctuations. 



counts per channel, thereby allowing x 1 statistics to be used. 
The data from the three cameras were modelled simultaneously 
using the XSPEC code, version 11.3.0. Spectral fitting is per- 
formed in the [0.5-8] keV band. The spectra are modelled using 
a simple, single-temperature model (MEKAL plasma emission 
code in XSPEC) with the absorbing column density fixed to the 
Galactic value (JV H = 4.85 X 10 20 cm~ 2 , Dickey & Lockman 
1990). The free parameters in this model are the temperature 
kT, metallicity Z (measured relative to the solar values) and 
normalization (emission measure). 

The best-fitting parameter values and 90% confidence lev- 
els derived from the fits to the annular spectra are summarized 
in Table 1. The projected temperature profile determined with 
this model is shown in Fig. [5] The temperature rises from a 
mean value of 8.8 + 0.3 keV within 115 kpc to kT = 11.1 ± 0.4 
keV over the 0. 1-0.5 Mpc region, then declines down to a mean 
value of 6.0^ 2 g keV in the outer regions (1.0- 1.7 Mpc). InFig.|5] 
we also show for comparison the projected temperature profile 
measured by Chandra (Allen et al. 2002). We note that while 
the general trend observed by the two satellites is consistent, 
there are some discrepancies in the measurements of the ab- 
solute temperature values. The discrepancy between Chandra 
and XMM temperature profile has been found in other clusters 
of galaxies (e.g. A1835, Schmidt et al. 2001, Majerowicz et al. 
2002), and can be partially due to the effect of the XMM PSF 
(see Markevitch 2002). A fit with the same model to the data 
for the SE quadrant between radii 50-200 kpc yields a best- 
fitting temperature kT = 13.3 +1.0 keV. In other directions, 
the mean value is kT — 11 .Of ke V. The metallicity profile 
derived with the single-temperature model is consistent with 
being constant, with an overall value of Z = 0.26 + 0.04Zq. 
However, as shown in Table 2, the structure in the innermost 
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Fig. 5. Circles: the projected X-ray gas temperature profile 
(and errors at 90% confidence levels) measured from XMM 
data in the [0.5-8] keV energy band. Triangles: the projected 
X-ray gas temperature profile (and la errors) measured from 
Chandra data in the [0.5-7] keV energy band (Allen et al. 
2002). The data of the SE quadrant are excluded in both pro- 
files. 

bin is better modelled by multi-temperature models having 
higher metallicities. 

Within the radius of 1.7 Mpc (~ 5'), a fit to the full 360° 
data gives an overall kT = 10.0 + 0.3 keV, Z = 0.26 + 0.03Z© 
and L x (2-10 keV) = 6.0 + 0.1 x 10 45 erg s -1 . These values are 
in agreement with ROSAT and ASCA results (Schindler et al. 
1997). 

6. Mass determination 

The total gravitating mass distribution shown in Fig. |6] (solid 
line) was calculated under the usual assumptions of hydro- 
static equilibrium and spherical symmetry using the depro- 
jected density distribution calculated from the parameters of 
the y6-model derived in Sect. [3] Only data beyond 30" (~ 175 
kpc) are considered: in the central bins the temperature as esti- 
mated in Sect.|5]is affected by the XMM PSF and projection ef- 
fects, while for the outer regions these effects can be neglected 
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(e.g. Kaastra et al. 2004). Within 1 Mpc we find a total mass 
of 1.0 + 0.2 x 10 15 M Q , in agreement with Chandra (Allen et 
al. 2002) and weak lensing analysis (Fischer & Tyson 1997) 
results and slightly higher than that derived by ROSAT/ASCA 
(Schindler et al. 1997). In Fig.[6]we also show for comparison 
(dashed line) the mass profile derived by assuming a constant 
temperature of 9.5 keV. 




radius (arcmin) 

Fig. 6. Solid line: Profile of the integrated total mass. Dashed 
line: Profile of the integrated total mass calculated assuming a 
constant temperature of 9.5 keV. Dotted line: Error on the mass 
calculation coming from the temperature measurement. 



7. Cooling core analysis 

We accumulate the spectrum in the central 30" (~175 kpc) by 
excluding the data for the SE quadrant. We use three different 
spectral models. Model A is the MEKAL model already used in 
Sect. [5] Model B includes a single temperature component plus 
an isobaric multi-phase component (MEKAL + MKCFLOW 
in XSPEC), where the minimum temperature, kT\ ovi , and the 
normalization of the multi-phase component, Normi ow = M, 
are additional free parameters. Finally, in model C the con- 
stant pressure cooling flow is replaced by a second isother- 
mal emission component (MEKAL + MEKAL in XSPEC). As 
for model B, this model has 2 additional free parameters with 
respect to model A: the temperature, kT\ OVJ , and the normal- 
ization, Normi ow , of the second component. The results, sum- 
marized in Table 2, show that the statistical improvements ob- 
tained by introducing an additional emission component (mod- 
els B or C) compared to the single-temperature model (model 
A) are significant at more than the 99% level according to the 
F-test, although the temperature of the hot gas is unrealisti- 
cally high. With our data, however, we cannot distinguish be- 
tween the two multi-phase models. This means that the extra 
emission component can be equally well modelled either as a 
cooling flow or a second isothermal emission component. We 
note that the fit with the cooling flow model sets tight con- 
straints on the existence of a minimum temperature (~ 2 keV). 
The nominal mass deposition rate in this empirical model is 



Table 2. The best-fit parameter values and 90% confi- 
dence limits from the spectral analysis in the central 0-30 
region. Temperatures (kT) are in keV, metallicities (Z) as 
a fraction of the solar value and normalizations in units 
of lO- u n e n p V/4nD A (l + z) 2 as done in XSPEC (for the 
MKCFLOW model the normalization is parameterized in terms 
of the mass deposition rate M, in Mq yr~'). 
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~ 1900 Mq yr _1 . This extremely high M, in agreement with 
Chandra results (Allen et al. 2002), is exceptional for distant 
cooling core clusters, and its implications on the evolution of 
cooling core clusters will be investigated in a forthcoming pa- 
per (Gitti et al. in prep.). 

8. Discussion and Conclusions 

The XMM-Newton observation of RX J1347.5-1 145 confirms 
that it is, with a luminosity L x = 6.0 + 0.1 x 10 45 erg s _1 (2-10 
keV energy band), the most X-ray-luminous cluster discovered 
to date. RX J1347.5-1 145 is a hot cluster (overall temperature: 
kT = 10.0 + 0.3 keV), not isothermal: the temperature profile 
shows the presence of a cool core and a decline of the tem- 
perature in the outer regions. The temperature map reveals a 
complex structure and identifies a relatively hot region at radii 
50-200 kpc to the SE of the main X-ray peak. This hot region 
is found at the same position as the subclump seen in the X- 
ray image. The higher temperature and enhanced X-ray surface 
brightness in the SE quadrant indicate that there is probably an 
ongoing subcluster merger event. On the other hand, excluding 
the data of the SE quadrant the cluster appears relatively re- 
laxed and the presence of a cooling core, with an exceptional 
high value of the mass accretion rate (M ~ 1900 Mq yr _1 ) in- 
dicates that the cluster could have evolved without any distur- 
bances for a relatively long interval. Therefore, the dynamical 
state of RX J1347.5— 1 145 appears very complex and will be 
presented in detail in a forthcoming paper (Gitti et al. in prep.). 
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